Bound pair states beyond the condensate for Fermi systems below T c : 
the pseudogap as a necessary condition 
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As is known, the 1/q 2 theorem of Bogohubov asserts that the mean density of the fermion pair 
states with the total momentum obeys the inequality n q > C/q 2 (q — > 0) in the case of the Fermi 
system taken at nonzero temperature and in the superconducting state provided the interaction 
term of its Hamiltonian is locally gauge invariant. With the principle of correlation weakening it is 
proved in this paper that the reason for the mentioned singular behaviour of n q is the presence of 
the bound states of particle pairs with nonzero total momenta. Thus, below the temperature of the 
superconducting phase transition there always exist the bound states of the fermion couples beyond 
the pair condensate. If the pseudogap observed in the normal phase of the high-T c superconductors 
is stipulated by the presence of the electron bound pairs, then the derived result suggests, in a 
model-independent manner, that the pseudogap survives below T c . 



PACS numbers: 74.20.-z, 05.30.Fk, 05.30.-d 



I. INTRODUCTION 



At present the pseudogap is well established to be in 
the spectrum of the elementary excitations of undoped 
and optimally doped-, high-T c superconductors (for ex- 
ample see the reviewu). The presence of the pseudogap 
implies that the electron subsystem in the normal phase 
is not the Fermi liquid and, so, theoretical explanation 
of the pseudogap is recognized as the key point of under- 
standing the phenomenon of the high-T c superconduc- 
tivityB'0. There are a great number of various theoretical 
approaches of investigating this problem. Two of them 
considered below are especially interesting in the context 
of this paper. 

The pseudogap can be associated with the presence of 
the local pairing correlations without phase coherence. 
The idea of this approach assuming the singlet pairing of 
fermions without the phase coherence, as applied to t 
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high-r c superconductivity, has been proposed in 
The more radical model of Alexandrov and Mot 
erates with, say, preformed bosons (bipolarons) existing 
in the system above T Cl the pseudogap being treated as 
coming from the binding energy of a bipolaron (of the 
order of a few hundred K). This model dates back to 
the Schafroth's ideas according to which the supercon- 
ductivity is a result of the Bose-Einstcin condensation of 
the bound pairs of electrons localized in the space and 
appearing in the system before the condensational. 

The concept of a bound state of two particles in 
medium can consistently be formulated with the reduced 
density matrix of the second order (2-matrix)t3. Indeed, 
the system of two particles is a subsystem of that of N 
particles. So, its state is not pure even in the situa- 
tion when all the system of interest has a wave function. 
In general a subsystem is specified by the density ma- 



trix (see, e.g. Ref.ETj. In particular, the reduced density 
matrix of the second order is of use when a noncoher- 
ent superposition of the pure states of two particles is 
relevant rather than any wave function. If among these 
states there exist bound ones, then a part of_particles 
of the system involved form bound pair statesEI. In the 
superconducting phase a macroscopical number of parti- 
cle pairs Nq occupy the same bound state, i.e. there is 
the condensate of pairs at which the ratio Nq/V = no 
is constant in the thermodynamic limit V — > oo. In the 
space-uniform case the condensate is formed by the pairs 
with the zero total momentum hq = 0, the binding en- 
ergy Eb of these-pairs being just the value of the super- 
conducting gaplilJ. The bound particle pairs beyond the 
condensate are characterized by the continuous distribu- 
tion over the total momentum of a couplalJ . The couple 
like these must also have the finite binding energy £fe(q) 
that, due to the continuity argument, should tend to e& 
when q — ► 0. If these bound particle pairs are 'hard' 
clusters, like in the theory of Alexandrov and Mott, then 
one may consider that the quantity £fc(q) is practically 
independent of q. The binding energy £&(q) is just the 
pseudogap, which manifests itself in the normal phase 
when the bound couples survive at T > T c . 

In the BCS-theory there are gp-.bound pair states be- 
yond the condensate absolutelyQii 2 ] (see below) which is 
a consequence of the violatkm of the local gauge invari- 
ance (see, for example, Ref.li3). 

In this paper we shall prove in a model-independent 
manner that the existence of the condensate of the bound 
pair states (BCS-pairs) implies the presence of the bound 
couples beyond the condensate (Schafroth's pairs). Em- 
phasize, that we do not specify the size of the pairs. If 
it is much more than the mean distance between par- 
ticles (the condensate pairs in the BCS-model), then, 
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following BogoliubovB, one may call these pairs 'quasi- 
molecules'. If the radius of the bound particle couples is 
of the order of the mean distance between particles or, 
even, less (the Schafroth-Alexandrov-Mott approach), 
then one may speak about an ordinary molecules. The 
proof is based on the well-known 1/q 2 theorem of Bogoli- 
ubov for the Fermi systemQ which is valid in the space- 
uniform case and under the condition of the local gauge 
invariance of the interaction term of the system Hamil- 
tonian. 

The present article is organized as follows. In section |0] 
the concept of in-medium wave functions of fermion pairs 
is considered. The properties of the pair condensate are 
discussed in the third section. At last, the proof concern- 
ing the noncondensed bound pairs of fcrmions is given in 
section |^ of the paper. 



II. THE CONCEPT OF PAIR WAVE FUNCTIONS 
FOR FERMIONS 

Thus, let us consider a homogeneous Fermi system of 
N particles with the spin s = 1/2 at nonzero tempera- 
tures. Suppose that the total momentum and spin of the 
system are conserved quantities. Let the forces exerted 
by fermions on each other be described with the two- 
particle interaction potential depending on the relative 
distance between them and, may be, on the spin vari- 
ables like in the case of various effective Hamiltonians. 
A state of the whole system is specified by the density 
matrix corresponding to the canonical Gibbs ensemble: 



p = exp 



H 



Tr exp I — 



H 



(1) 



where H is the system Hamiltoniantfl. In this case thp 
2-matrix is represented in the form (see, for example,£3) 

= N ,*_ 1} {^{x l )i,\x 2 )^{x' 2 )i,{x' 1 )), (2) 

where (• • •) = Tr (• • • jo) stands for the average over the 
state (|l[); x = (r, a) represents the space coordinates r 
and spin z-projection a = ±1/2; ip^(x), ip(x) are the 
field Fermi operators. The 2-matrix obeys the normal- 
ization condition 



J dx\dx 2 p2{xi,X2',Xl,X2) = 1, 



(3) 



here J dx ■ ■ ■ = y] J d 3 r ■ ■ ■ and integration is fulfilled 
over the volume V. Therefore, the 2-matrix (||) has the 
asymptotic behaviour 1/V 2 when V — + oo, n — N/V — 
const. So, it is more convenient to deal with the pair 
correlation function F 2 differing by a norm from p 2 : 

F 2 (x 1 ,x 2 ;x' 1 ,x' 2 ) = (^(xi)^(x 2 )i>(x' 2 )iP(x[)). (4) 



The boundary conditions for ^0 follow from the prin- 
ciple of |the correlation weakening at macroscopical sep- 
arationsQ: 

(^1)^2)^2)^1)) - 

(^( Xl )^(x 2 )) (v>(4)V>04)) (5) 

when 

ri — r 2 = const, — r' 2 = const, \r[ — ri| ^00; (6) 

(■0 1 '{Xl)ll) ] ' (X 2 )lp(x' 2 )lp(x' 1 )) -> 

(V^iMz'i)) {^{x 2 )^{x' 2 )) (7) 

when 

ri — r[ = const, r 2 — r' 2 = const, |ri — r 2 | — > 00. (8) 

As the kernel (^) is a non-negative Hermitian operator 
acting on the two-particle wave functions ^(xi,x 2 ), we 
can expand it in the orthonormal set of its eigcnfunctions 
(EF): 

F 2 {x 1 ,x 2 ;x' 1 ,x' 2 ) =22n i/ %I)*(xx,x 2 )'iIj v (x' 1 ,x 2 ), (9) 

v 

where 



dx\dx 2 \^ v {xi^x 2 )X 



(10) 



Eigenfunctions ipv(xi,X2), which at the same time are 
EF of 2-matrix (§), are called the pair wave functions, 
or PWF. 

With (|j), (||) and ( ^0|) one can be convinced that 



dxidx 2 F 2 (xi, x 2 ; xi, x 2 ) = (N 2 - N) 
= N(N- 1) 



Therefore, the non-negative quantity N v can be inter- 
preted as the mean number of the pairs in the state v, any 
pair being doubly taken. The ratio w v = N„/{N(N — 1)} 
is the probability of observing a particle pair in the pure 
state with the wave function ip v (xi, x 2 ). Here, as one 
might expect, ^2 v w v = 1. 

It follows from the definition (^) that 

F 2 (xi,x 2 ;x' 1 ,x' 2 ) = -F2(xi,X2',x , 2,x , 1 ) = 
= -F 2 (x 2 ,xx;x' 1 ,x 2 ). 

So, tp v (xi, x 2 ) = —ip u (x 2 ,xi), i.e. PWF for fermions, as 
usual, are antisymmetric with respect to permutations of 
particles. 

In an equilibrium state the total pair momentum ?iq is 
a good quantum number for PWF provided that the total 
momentum of thej-whole system is a conserving quantity 
(see proof in Ref.ll3). The same is correct for the total 
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spin S of a particle pair if there is no magnetic ordering!* 
So, the index v can be represented as v = (w, q, S), where 
u! stands for other quantum numbers. As to the PWF, 
they can be written as 

cxp{iq(ri +r 2 )/2} 
= ipu,q,s(ri - r 2 ,cri, cr 2 ) j= . (11) 



Then expression (^J) has the form 



V 



F2(xi,X2;x' 1 ,x' 2 ) = ^2 NuJ y' S - C, q ,g( r i - r2,gi,o- 2 ) 



cj,q,5 



xV>o;,q,s(ri-r2, (Ti,(T2) expji^ri+r'a-ri-ra) I . (12) 

For the wave function "0c*j,q.s( r : cr i; cr 2) which can be 
interpreted as the wave function of a particle pair in the 
center-of-mass system, from ( |lO| ) and (|ll]) we obtain 



d 3 ri^q,s(r,eri,(T 2 )| 2 = 1. 



(13) 



It can be related to either discrete or continuous spectra. 
In the former case 



VVq,s(r, CTi, (T 2 ) -> 



(14) 



when r — > oo and, so, we deal with the sector of bound 
states of particle pairs. The latter variant implies 

VVq,s(r,eri,cr 2 ) -> x s , ms (ffi, <t 2 )v / 2 | ^j^j (15) 

for r — * oo. This is a 'dissociated', or scattering, pair 
state corresponding to the relative motion with the mo- 
mentum hp. Here XS,m s i. cr i-> cr 2) is the spin part of the 
pair wave function (spinor), mj being the z -projection of 
the total pair spin S. When S = 0, ms = (the singlet 
state) one should take cos(pr). For 5=1, ms = —1, 0, 1 
(the triplet state) one should use sin(pr). Remark that 
in the situation when the fermion interaction does not 
depend on spin variables, the spin and space parts of the 
PWF can be separated from one another not only when 
r — > oo but also for any r. 

In the case of ( |l4| ) to = i, where i stands for the dis- 
crete index enumerating the bound pair states. Let us 
denote ip u> q,s(r, <ti, <7 2 ) = <Pq,s,i(r,o-i,0- 2 ), so that 



E 



J d 3 r \(pi,q,s(r, eri,er 2 )| 2 = 1. 



(16) 



In the situation of ( |15[ ) u) = (p, ms)- Here it is convenient 

to introduce Vw,q,s( r > <ri> a i) — Vp,q, S,m s 
From (||) it follows that 

7/ J2 J^r\ip PA ,3,m s {r 1 a 1 ,<T 2 )\ 3 = l. (17) 



Now, with the variables 

R=(ri+r 2 )/2, r = ri - r 2 



(18) 



and, respectively R' and r', the expression (jlj) is rewrit- 
ten as 



F 2 {xi,x 2 ;x'^x' 2 ) 



E Nq,S,i 
v r,! 

1,S,i 



( Pq,S,i( T ' <Jl ' (72 ) 



><¥'q 1 s,i(r / .^i ? ^)exp{iq(R / - R)} 



x 



Jv p.q.S,m s * / \ 

N — ^ Vp,q,S,m s ( r . CT 1^2) 

p,q,S,ms 

V P ,q,s,m s (r',cri,(r 2 )exp{iq(R' - R)} . 



(19) 



In the thermodynamic limit all the summations over mo- 
menta can be replaced by the corresponding integrals: 

F 2 (x 1 ,X 2 ;x' 1 ,x' 2 ) = ^2 / <P9 w S,i(<l) ^q,S,i( r )°'lj a '2) 
S,i J 

xp q ,s,«(r', (ri, o- 2 ) exp {iq(R' - R)} 

+ E / d3 P d3 9 W -S,m s (P,q) ¥'p,q,S',TO S ( r ) 'l ! cr 2) 

c ^ 



:yp,q,S,m s (r',cr 1 ,f7 2 )exp{iq(R' - R)} 



(20) 



Thus, from equations (|19|) and ( ]20| ) we can see that 
V ws,i(q)d 3 q is the number of the bound particle pairs 
with the spin S, in the state i and with the total couple 
momentum frq located in the infinitesimal volume d 3 q. 
Respectively, V 2 ws,m s (p, q)d 3 pd 3 q stands for the num- 
ber of the 'dissociated' particle pairs in the state (S, ms) 
with the relative momentum hp and total momentum h(\ 
located in the infinitely small volumes d 3 p and d 3 q. 

In the center-of-mass system the replacement p — > 
— p, <j\ — » er 2 , <T2 — > (7i corresponds to the permutation 
of particles. So, the following symmetric relations take 
place: 



ws, ms (p, q) = w s , ms (-P) q)> 



(21) 



^p,q,5,ms (l") (7l 5 (7 2 ) = - ^p,q,S,m s (~ r , CT 2 , (7 1 ) 

= -<^-p,q,S,m S (r, (72,(71). (22) 

As an example, let us consider the expansion of F 2 in 
terms of PWF for the BCS-model. Taken with an accu- 
racy to the asymptotically small quantities, the Hamilto- 
nian in the BCS-approach is represented as the quadratic 
form of the Fermi operatorstJ that can be diagonalized 
with the Bogoliubov transformation. Therefore, one is 
abj£j to use the theorem of Wick, Bloch and Dc Domini- 

F 2 (xi,x 2 ;x[,x 2 ) = (V' t (xi)ip^ (x 2 )i/j(x' 2 )ip(x' 1 )) 
= (^( Xl )^(x 2 )) (V(a£Mzi)> + (^(ziMzi)) 
x(^(x 2 )i>(x' 2 )) - {^{x x )^{x' 2 )) (^(x^x'J). (23) 
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Further, for the 'normal' averages we have 
(^(anMsi)) = (^('i.o-i^W.Oi)) 
= / ^ynWexp^kK-rOlA^-ai), (24) 

where n(k) = (a^ ff ak iff ) gives the distribution of 
fcrmions over momenta; and we introduced the function 



A(a) 



0, (7^0, 

1, (7 = 0. 



'Anomalous' averages are given by 

(^( a;i )^K)) = (^(r 1 ,a 1 )^(r' 1 ,ai)) = 
f d 3 k 

=y7^T3( a k^i a -k,-o-i) ex P{* k ( r l - ri)}A((Ti+cri). (25) 

In the BCS-model, the quantity (ak,<r&-k,-er) can be rep- 
resented in the following form 



(flk^a-k-ff) = V^o#) 



sign(cr) 
V2 ' 



with y(fc) obeying the normalization condition 
' /!/ '' l^)| 2 = l. 



(26) 



(27) 



Remark that one can consider ip(k) as a real quantity be- 
cause it can be made real with the corresponding phase 
transformation of the operators and a^. Now, Eqs. 



(|2J), (|2J) and (J26J) allow us to rewrite in the follow- 
ing form: 

F 2 (x 1 ,x 2 ;x' 1 ,x' 2 ) = n ip(r)xo,o{(Ti, <T 2 )(p(r')xo,o(v'i, v' 2 ) 
+ ^ l ~d 3 pd 3 q n f(i 

S,ms 



(2tt) 6 "V2 



c<Pp,s(r')xs,m S (o-i,(72)exp{iq(R' - R)} . 



(28) 



Here <p(r) is the Fourier transform of ip(k), for y P) s(r) 
we have 



V P ,s( r ) 



V^cos(pr), 5 = 0, 



v / 2sin(pr), 5=1. 
Respectively, the spinor x stands for 

XS,m s (ci, cr 2 ) 

A((7i + cr 2 )sign(c7i)/V2, 5 = 0, m s = 0; 

e(-ffi)0(-(r 2 ), 5=1, m s = -l; 

A((7i + «7 2 )/V2, 5=1, m s = 0; 

G(<7i)e(a 2 ), 5=1, m s = l. 

Here 



(29) 



(30) 



6(a) 



1, a > 0, 
0, cr < 0. 



With (|27j), ( p9| ) and j30| ) one can easily be convinced 
that the normalization relations (jl^) and ( |l7| ) are satis- 
fied. Within the BCS-model w s ,i{q) = A(5)A(i) n <$(q) 
(<5(q)) is the ^-function), i.e. all the bound particle pairs 
are condensed. 



III. PROPERTIES OF THE CONDENSATE OF 
PAIRS 



Let us demonstrate in the most general case that if the 
'anomalous' average (^(x^tp^x^) is not equal to zero (off 
diagonal long-range order) then the distribution function 
Ws,i(q) acquires the 5— functional singularity correspond- 
ing to some indices Sq and iq or, in other words, the ratio 
^Vq,s ,i /F in the first sum of ( |l^ ) does not vanish in the 
thermodynamic limit: 

wsAl) = n <5(q)A(5 - 5 )A(i - i„) + ufe.^q), (31) 



where ws.i(q) is the regular part of (J3lJ) giving the 
bound-pair distribution over nonzero momenta. 

To do this, let us take the limit relation (JsJ) and rewrite 
it with the variables ( |l8| ) in the form 

F 2 (xx,x 2 \x\,x' 2 ) -> (4>\xi)ip\x 2 )) 

x(ip(x' 2 )ip(x[)) =n ip*(r,a 1 ,a 2 )v(-r',a' 1 ,a 2 ), (32) 

where the functions ip*{v, oi, <7 2 ) and y(r', a' 2 ) are in- 
troduced in such a way that the normalization condition 
( |l6| ) should be fulfilled. This can always be done because 
according to the principle of correlation weakeningu 

(V>(ziM*2)) -> M*i)> (#* 2 )> = o 

when r — > oo (see Ref.0). Expression ( |32] ) is exactly the 
contribution of the first singular term o f (fn] ) into (|2Cj). 
The contribution of the regular part of (|3l|) and that of 
the 'dissociated' pair states into (^0|) are infinitely small 
in the situation of (||) according to the Riemann's theo- 
rems because 



ri +r 2 



IR' - RI 



Remark that the pair distribution over the scattering 
states ( f2l| ) does not contain 8— functional terms. Indeed, 
in the opposite case they would lead to the condensate 
of thepane-particle states like in the situation of the Bose 
liquidEZl which is impossible for the Fermi systems. 

Eq. (H|) allows us to treat the 'anomalous' averages 
as the wave functions of the condensed pairs of fermions, 
(of course, with an accuracy to the normalizing factor). 
For the density of the pairs like these Eq.([l6|) and ( |3lq ) 
gives 



4 



n 



V 



= £ / ^|<V>(r,oi)#},o 2 )>| 2 = 



E 



(27T) 



(33) 



where it has been taken into account that the total mo- 
mentum of the system and ^-component of its spin are 
conserved quantities. Keeping in mind these integrals of 
the motion, one could expect that in the most general 
case the wave function of the condensed pairs should be 
written as 

(p{r,ai, cr 2 ) = -^(ip(r 1 ,a 1 )tfj(r 2 ,a 2 )} 



d 3 k A((Ti + a 2 ) 
3 (fflk.o-ia-k.-crj f= exp(zkr) 



(2tt) 

1 A(o-i + a 2 ) 
V2 



(<p s (r)sign(oi) + ^t(r)), 



(34) 



where (f s (r) = ip s (—r) and <£>t(r) = — ft{— r). Accord- 
ing to Eq. (|3(]) the first term in ( |34| ) corresponds to the 
singlet and the second, to the triplet components of the 
wave function of the condensed fermions. However, (|34|) 
is not quite correct because the total pair spin should be 
an integral of the motion, even in the situation with the 
spin-dependent interaction between fermions. Therefore, 
we are not able to obtain a superposition of the singlet 
and triplet states. Instead, in ([54]) one should select ei- 
ther (p s (r) ^ 0, (p t (r) = or ip s (r) = 0, <p t (r) ^ 0. So, 
we have: 

f fs(r)xo,o{<T 1 ,a 2 )/y/^o, 
<p{T, o-i, cr 2 ) = <^ (35) 
I <Pt(r)Xi,o(o-i,o 2 )/v / n^. 



The phase coherence takes place for the condensed 
bound pairs due to the uncertainty relation AipANg ~ 1 
for the phase tp and number of the bound fermion pairs 
No = Aq = o,So,«o m the state (q = 0, So,io)- In the 
thermodynamic limit the macroscopical occupation of 
this state results in ANq oc — > oo and, therefore, 
A<p — ► 0. For the bound pair states beyond the conden- 
sate Nq^s,i is limited above even for V — » oo. Thus these 
states are not correlated with respect to the phase. 

Remark that the total number of the bound particle 
pairs (condensed and not) 



N b = V 



S,i 



d 3 qw s ,i(q) 



is proportional to the total number of particles N. In 
particular, there is the inequality for the number of the 
condensed bound pair statesc3 



N < N. 



(36) 



It should be emphasized that the inequality (|3(j) is not 
trivial. One can consider, for example, a dilute gas of to— 
particle molecules. In this case we have 2V& = (m — l)iV, 
thus, one can obtain Nb > N provided m > 3. 



In the space-uniform case we can readily find relation. 
d36| ) with the inequality of Cauchy-Schwarz-Bogoliubovtl 

\{AB)\ 2 < {AA} ){&B). 
Indeed, assuming A — a^.a and B — a_k,-er we arrive at 

|(ok,ff a-k.-cr>| 2 < (ak, CT a kiCT ) (al k) _ CT a-k,-<r) 

= (1 - n{k)) n(k). 

Then, from (|3^) we derive 
N 



n = 



V V 



a-k. 



2 2 N 

k k 

It is interesting to note that n(k) — n 2 (k) — 

((a k ,o- a k,o-) 2 ) - {a^a^a) 2 = D(n(k)) is nothing else but 
the mean square deviation of the occupation number of 
the (k, a) one-particle state. So, the stronger inequality 



no< v J2 D H k )) 



(37) 



demonstrates that the number of the condensed pairs is 
tightly connected with the 'wash-out' of the Fermi sur- 
face. In the BCS-model at zero temperature 



no_ 

n Ep 



< 1 



because the bound pairs are formed by the particles lo- 
cated near the Fermi surface only. In general, no is the 
most 'reliable' order parameter of the superconducting 
phase transition. 



IV. THE BOGOLIUBOV 1/Q 2 THEOREM AND 
BOUND PAIR STATES BEYOND THE 
CONDENSATE 

Let us now prove with the principle of the correlation 
weakening that the distribution of the particle pairs over 
the 'scattering' states ws, ms (p, q) is expressed in terms 
of the occupation numbers of the one-particle states 
n(k) = (a k g.ak.0-). Indeed, on the one hand, in the lim- 
iting situation of (||) we have the relation (^), which can 
be written as 



F 2 {xi,x 2 \x' 1 ,x' 2 ) 
d 3 P i 



(2tt) 3 
d 3 P2 



(27T) 



n{pi) exp(jpi(ri - n))A(oi - <j[) 
3 n(p 2 )exp(ip 2 (r 2 - r 2 ))A(o 2 - <j' 2 ) 



(27T) 



exp(ip(r' — r)) 



x exp(iq(R' - R))A(oi - a[)A((7 2 - a' 2 ), 



(38) 
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where, passing to the last equality, we introduced the 
new variables q = pi + p 2 and p = (pi — P2V2 and used 
notations (|l8|). On the other hand, when (||) is true, we 
have 



r 2 -ri 



00, r 



|r + r'| — > 00, R' R = const, r' — r = const. 

Therefore, it follows from |L4|), (|l|) and @ that in the 
limiting case (H) we have 

F 2 (x 1 ,x 2 ;x' 11 x' 2 ) -> 

^ y d 3 qd 3 p ius )ros (p,q)v?p,s(r)ipp,s(r') 

S,ms 

XXS,m s (o-i,cr2)XS,m s (o-i,cr2) ex P(*q( R ' ~ R )) ) ( 39 ) 



where we used notations (29) and (|30|). Further, the 
Riemann's thcoremE3 used while integrating over p and 
relation (21) allow us to rewrite (|39|) as 



F 2 {x 1 ,x 2 ;x' 1 ,x 2 ) -> 

d 3 (? cZ 3 p w s , ms (p, q)xs,m s (ci , ff 2 ) 

S,ms 

XXS,m S ((r[, o' 2 ) exp (ip(r' - r)) exp (iq(R' - R)) . (40) 



The right-hand side of Eq. (|38|) is equal to that of (|4G|) 
at all the values of the spin variables and space ones 
f = r' — r and R = R' R. Taking into account the 
completeness of the set of the spin functions (|30|) 



^2xS,m s (vi, V2)XS,ms(v'l, V 2 ) = A((7i - <t[) A((7 2 - <T 2 ) , 
S.m s 

we derive the following equality: 

n(q/2 + p) u(q/2 - p) 



ws, ms (p,q) 



(2tt) 



Thus, in the thermodynamic limit one can write 

Np^s, ms = «(q/ 2 + p) «(q/ 2 - p)- 



(41) 



(42) 



As it is seen, when there is no magnetic ordering (it is 
obviously true for the superconducting phase) , the func- 
tion of the pair distribution over the 'dissociated' states 
is independent of the quantum numbers S, ms- 

It is now easy to prove that the pair condensate must 
always be accompanied by the presence of the noncon- 
densed fermion pairs: u>s,i(q) 7^ in ( |3l] ) if n 7^ 0. Let 
the interaction energy of the system be invariant with re- 
spect to the local gauge transformation of the field Fermi 
operatorsE-3 



i/j(r, a) -> tp(r, a) exp (ix(r)) , 
V> f (r, a) -> V + (r, °) exp (-i%(r)) . 



In this case the 1 jq 2 theorem of Bogoliubov for the Fermi 
systemsQ is valid which asserts that in the presence of the 
pair condensate we have the inequality for sufficiently 
small q 

C 

uj,S M q 2 



where N^.q.s appears in Eq. ([12|) and u is the set of the 
quantum numbers corresponding to both the continuous 
spectrum (uj ~JjP, ms)) and the discrete one (uj — i)cB. 
However, Eq. (f42j) results in 



N, 



p,q, 5,771; 



< 1 



because n(k) < 1 for fermions. Therefore, we have the 
only possibility at which the singularity 1 / q 2 appears due 
to the noncondensed bound pairs. It is reasonable to ex- 
pect that these pairs have the quantum numbers of the 
condensate couples So,«o : 



C 

ws ,*o(q) > — • 

q z 



(44) 



(43) 



The BCS-model is not locally gauge invariant which 
results in absence of the noncondensed bound pairs: 
S's.i(q) — 0. It is important to note in this connection 
that the bound pair states beyond the condensate may 
play a noticeable role in calculating the gauge-invariant 
response of the system to the electromagnetic fields. 

We have proved that the noncondensed bound pairs 
coexist with the condensed ones at T < T c . So, any the- 
ory ignoring the noncondensed bound pairs of fermions 
is not fully consistent. Remark that the distribution of 
the bound fermion pairs over the center-of-mass mo- 
menta obeys the inequality ( ^i|) with C oc ksTno (see 
Rcf.El). The distribution of the particles over momenta 
in the Bose gas w(q) = n(q)/(2ir) 3 answers, at small q, 
the similar relation w(q) > C" /q 2 with C" oc fcgT' 
(here no denotes the density of the condensed bosons)! 
Therefore, there are fundamental parallels between the 
Bose gas and the considered subsystem of the fermion 
bound pairs. And these parallels are not only reduced to 
agreement between the fermion-pair statistics and the 
Bose one. Following this analogy, we can expect that 
the bound fermion pairs exist even at T > T c (appar- 
ently, in some temperature interval T c < T < T*, in 
spite of the disappearance of the 1 / g 2 -singularity) . Thus, 
it looks as if any superconducting phase transition is a 
particular case of the Bose-Einstein condensation. This 
conclusion can be of interest in the context of the dis- 
cussion concerning different approaches pf. investigating 
the high-T c superconductivity (see Refs.c3'E3). Remark 
that possible experimental consequences of the existence 
of fermion boixad pairs beyond the condensate can be 
found in papertj in the case of neutral Fermi systems. 

The space-uniform character of the Fermi system is of 
use in the proof given above. Electrons in the crystalline 
field, of course, can not be treated on the same level. 
However, for g -t (large wave lengths) a crystalline 



G 



lattice can be considered as continuum. Therefore, the 
derived result remains correct in this case. 

Emphasize that the bound pair states can fully be a 
result of the collective effects. Indeed, as it was demon- 
strated by CooperES, an arbitrary small attraction be- 
tween electrons leads to forming the condensate of the 
bound electron pairs. Hence, if we considered a suffi- 
ciently shallow well as the two-fermion interaction poten- 
tial, we would observe formation of the condensed and, 
according to the obtained result, noncondensed pairs at 
low temperatures. However, the well can be chosen in 
such a way as to prevent the bound states of two 'bare' 
fcrmions from appearing within the ordinary two-particle 
problem. 

At last, it is important to make one more remark on 
the connection between the Xiq 2 theorem of Bogoliubov 
and the Goldstone theoremCJ. As it has been demon- 
strated in RefE, existence of the Goldstone mode in the 
Bose system results from the Bogoliubov theorem pro- 
vided the mass operator k) is regular in the vicinity 
of the point lu = 0, k = 0. Let us emphasize that there 
are situations when the Bogoliubov theorem is valid while 
it is not the case for the Goldstone one. For example, in 
the case of neutral weakly interacting Bose gas the con- 
dition mentioned above for the mass operator is correct, 
and the Goldstone mode exists. On the contrary, for 
the charged Bose gas the mass operator is not regular 
at k = 0, and, thus, there is no Goldstone mode. The 
similar situation is realized for the Fermi systems (see, 
e.g. RefB). 



V. CONCLUSION 



Concluding, let us take notice of the main results once 
more. The reduced density matrix of the second order is a 
fundamental characteristic of a many-particle system, its 
eigenfunctions being the pure states of two particles se- 
lected in an arbitrary way. Appearance of the condensate 
of the bound pair states ( |33| ) implies the occurrence of the 
S— functional term in the distribution of the bound pairs 
over the momentum of the couple center of mass q (see 
Eq. (|3l|)). Using the space homogeneity of the system 
and the local gauge invariance (43) of the fermion inter- 
action, we have proved that there is the l/<7 2 -singularity 
in the distribution function u>s,i(q) provided that uq =/= 0. 
Thus, we refined the 1/q 2 theorem of Bogoliubov, having 
proved the singularity to appear in ws t i{q). Therefore, 
presence of the noncondensed bound pairs below T c is 
the necessary condition of superconductivity. 

A new simple proof of the Yang inequality for the Fermi 
systems (36) and its stronger variant (^) have also been 
derived as results of secondary importance. 

This work was supported by the RFBR Grant No. 97- 
02-16705. Discussions with V. V. Kabanov and V. B. 
Priezzhev are gratefully acknowledged. 
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